BL Lacertae was the target of an extensive multiwavelength monitoring campaign in the second half of 2000. In this paper, we are using leptonic and hadronic jet models to fit the observed broadband spectra and spectral variability patterns. We start out with global spectral models. Subsequently, we investigate various flaring scenarios for comparison with the observed short-term variability. For our leptonic jet model, we find that the short-term variability, in particular the optical and X-ray spectral variability, can be best represented with a flaring scenario dominated by a spectral-index change of the spectrum of ultrarelativistic electrons injected into the jet. Based on this result, a detailed model simulation, reproducing the observed optical and X-ray spectral variability and broadband SED of BL Lacertae simultaneously, is presented. Our leptonic modeling results are compared to fits using the hadronic synchrotron-proton blazar (SPB) model. That model can reproduce the observed SEDs of BL Lacertae in a scenario with muon-synchrotron dominated high-energy emission. It requires a significantly higher magnetic field than the leptonic model (~ 40 G vs. ~ 2 G in the leptonic model) and a lower Doppler factor associated with the bulk motion of the emission region (D ~ 8 vs. D ~ 18 in the leptonic model). The hadronic model predicts a significantly larger > 100 GeV flux than the leptonic models, well within the anticipated capabilities of VERITAS and MAGIC.